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I propound a non-linear generalization of the scalar-field Poisson equation of the form 
[(ip t z ifi_ i ) D ^ 2 ^ 1 ip k ]. k oc p, in curved D dimensional space. It is derivable from the Lagrangian den- 
sity L D — L° — Aptp, with Lj oc — (io, 1 co, i )" D/ ' 2 , and p the distribution of sources. Specializing 
to Euclidean spaces, where the field equation is V ■ (\X7tp\ D ~ 2 X7tp) oc p, I find that Lf is the only 
conformally invariant ^CI) Lagrangian in D dimensions, containing only first derivatives of <p, beside 
the free Lagrangian (Vco) 2 , which underlies the Laplace equation. When ip is coupled to the sources 
in the above manner, L D is left as the only CI Lagrangian. The symmetry is one's only recourse 
in solving this non-linear theory for some nontrivial configurations. Systems comprising N point 
charges are special and afford further application of the symmetry. In spite of the CI, the energy 
function for such a system is not invariant under conformal transformations of the charges positions. 
The anomalous transformation properties of the energy stem from effects of the self energies of the 
charges. It follows from these that the forces Fj on the charges q t at positions r* must satisfy certain 
constraints beside the vanishing of the net force and net moment: For example y^. -Fj must equal 
some given function of the charges. The constraints total (D + 1)(D + 2)/2, which tallies with the 
dimension of the conformal group in D dimensions. Among other things I use all these to derive 
exact expressions for the following quantities: 1. The general two-point-charge force. 2. The full 
potential field for two opposite charges ±g. 3. The energy function and the forces in any three-body 
configuration with zero total charge. 4. The few-body force for some special configurations. 5. The 
virial theorem for an arbitrary, bound, many-particle system relating the time-average kinetic energy 
to the particle charges. I also discuss briefly multi-scalar theories, theories with higher derivatives, 
and vector- and higher-form-potential theories. 

PACS numbers: 41.20.Cv, 03.50-z, 11.25.H 



I. INTRODUCTION 



It is a well-known and well-used fact that the Poisson equation, Ac/? oc p, for the potential ip produced by sources p, 
describes a conformally invariant (CI) theory in two dimensions: It is invariant under the angle-preserving coordinate 
transformations. In all dimensions it is linear in the field ip, and thus describes a "free" (non-self-interacting) field. 
Many of the special features of the D — 2 theory stem from its linearity, but many are underpinned by the conformal 
invariance. The Poisson equation in D > 2 dimensions is not CI (while the Laplace equation is-see section IIB). 

The Poisson equation describes many physical problems in linear media such as electrostatics, magnetostatics, 
steady-state diffusion and other potential flows in the presence of sources and sinks, and, of course, Newtonian 
gravity. It can be generalized to 

V • [m(|V</j|)V^] oc p, (1.1) 

to describe, for example, non-linear media with a response coefficient (dielectric constant, permeability, diffusion 
coefficient, etc.) that is a function of the field strength. An equation of this type, with different forms of p:(x), has 
been studied in different contexts. For example, as an effective-action approximation to Abelianized QCD Q, as a 
modification of Newtonian gravity to replace the dark-matter hypothesis for galactic systems ^ , and in the context 
of non- linear co mp osite media []4j. Some of the general properties of such theories are summarized and extended in 
p| . Equation ( p_.l| ) might serve as a model for many other nonlinear phenomena, such as electrodynamics in very 
strong fields. 

Here I point out that with the special choice of p(x) oc x D ~ 2 the theory is a natural generalization of the (linear) 
two-dimensional Poisson theory. The resulting non-linear theory is unique in certain regards. Foremost is its conformal 
invariance. This enables one to say much about the theory and its solutions-much beyond what is possible for the 
general case. The theory seems to be the only one derivable from a CI action that contains only first derivatives of the 
potential, with the source p coupled directly to the potential, i.e. with an interaction Lagrangian of the form pf(ip). 

In the modified dynamics discussed as an alternative to dark matter, phenomenology requires just this CI behavior 
in three dimensions in the limit of very small |Vcp| (see [£| ||). Our results here apply then in the large- distance limit 
of this theory. 
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In material media, non-linearities of the response coefficient appear at high values of \Vp\. Our results might then 
apply in the short-distance limit. So, for example, our results for point charges will be valid when charges are very 
near each other, and those for the fields at short distances from the sources. 

The present theory constitutes an instance of a highly non-linear theory that can be solved for non-trivial configu- 
rations due to the symmetry. 

I shall present two types of results: One concerns solutions for the potential field for various charge distributions 
obtained by conformal transformations from highly symmetric ones; this I do in section III, after discussing some 
general properties of the theory in section II. 

The other type of results concern systems of point charges. The dynamics of these is governed by an energy 
function that depends on the charges and their positions. It turns out that while the theory is invariant under 
conformal transformations, the energy-surprisingly perhaps-is not invariant under a conformal transformation of the 
positions of the point charges. This can be seen already in the two-dimensional case, which is exactly solvable, where 
the energy of a system of charges q t at positions r*, is E = J2i^j QiQjl n \ r i ~ r j\- Under a dilatation r s ; — > Ar^ we 
have E — > E + InX iiQj = E + (l/2)inA[(^ i qi) a — Yli if]} with a = 2. In the D-dimensional, non-linear case 
we do not, in general, have a closed expression for the energy. Still, we shall see that the energy transforms under 
dilatations as E — » E + KlnX, with K a function of the charges of the same form, with a value of the power a that 
depends on D. The non-trivial term in the transformation law comes from the behavior of the self energies of the 
charges under dilatation, including the fact that the self energy of a charge is logarithmic in its size scale. All this is 
rather transparent in the linear two-dimensional case. There is also an appropriate transformation law of the energy 
under inversions-the other conformal transformations (and of course, the energy is invariant under translations and 
rotations, which do not affect the self energies of charges). I discuss all this in section IV. Some of the applications 
to calculating energies and forces are brought in section V. In section VI, I discuss other field actions for ip, and 
demonstrate the uniqueness of L D as a CI Lagrangian. In section VII, I discuss multi-potential theories. In section 
VIII I generalize briefly to non-linear, CI extensions of Maxwellian electrodynamics in D > 4 dimensions. In the last 
section I bring brief comments on possible connections with quantum field theory. 



II. GENERAL PROPERTIES 

Via the equation 

V-{[(V^) 2 ] D / 2 - 1 V^} = a D Gp (2.1) 

a charge distribution p(r) in D-dimensional, Euclidean space gives rise to a potential <p(r). This field equation is 
derivable from the action 

S° = S? + S?^-J d°r -^-^ J [(V^) 2 ]^ 2 d°r. (2.2) 

Here, S° is the field action, S? is the interaction action, G is a coupling constant, and a D = 2(tt) d/2 /T(D/2) is the 
D-dimensional solid angle, introduced here for convenience. For G > 0, like charges attract, as in gravity; for G < 
they repel each other, as in electrostatics. The field equation has a unique solution in a volume V bounded by £ 
when either ip or the normal component of [(V '<p) 2 ] c '/ 2_1 V 'ip are dictated on £ (see e.g. M Jll]]). 

Two integral relations satisfied by solutions of the field equation were derived in |Tl| ] . The first applies for our class 
of theories when the total charge vanishes; it then tells us that 

^ / [(V^) 2 ]^ 2 d°r = -J pp d D r. (2.3) 

The second relation is an explicit expression for the virial integral V: 

V= / pr-Vp d D r= {dG)- 1 \GQ\ d = V D {Q) (2.4) 



[d = D/(D — 1)], which follows by writing V as a surface integral at infinity. The virial-which is shown below to 
control the response of the configuration's energy to rescaling-can then be written in terms of only the total charge 
of the system. 

While I shall work in Euclidean space with its specific conformal transformations, it is useful to formulate the 
problem for curved space. The covariant form of the action is 
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S° = I 9 1/2 PV d D r - ^J-^ J g 1 '*^,^) ' 2 d°r, (2.5) 
giving rise to the field equation 

l(g*^,^, J ) D/2 - 1 g km <P, k U = a D Gp. (2.6) 

Above, g i} is the metric, g' 3 its inverse, g = \det(g i:j )\, and summation over repeated indices is understood everywhere. 

The density, p, is defined so as to be a coordinate scalar: the charge within a volume V is S v g x ^ 2 p d D r. So, for 
example, for a system of point charges tfc at Ti 

p(r)=. 9 - 1 / 2 ^ (Z!( 5-(r-r 4 ). (2.7) 



Using usual derivatives instead of covariant ones eq. (|2.6D reads 

g- 1/2 [g 1/2 (g ij ^, } ) D/2 - l g km v, k U = a D G P . (2.8) 

The field stress tensor is defined as the functional derivative of the field action with respect to the metric; i.e., under 
a small change Sg i:j in the metric 

SSf = \j 'g 1/2 Sg^ d°r, (2.9) 

giving 

r« = -j^iv, V.0 D/ V - d^L), (2.io) 

which has a vanishing trace. (The metric does not appear in the interaction part-because g x ^ 2 p depends only on 
the charges degrees of freedom-which, thus, does not contribute to V.) The tracelessness results from the conformal 
invariance of the action Sf , as is well known (see below). For the Euclidean case the stress tensor becomes 

V = -(DaoGy^Vip^il- L>n®n), (2.11) 

with n = V(/?/|V</?|. 

In this flat case, the stress tensor gives the force on any volume V, bounded by the surface S on which p = as 



pVip d D r= V ■ ds. (2.12) 

J s 

(Compare with the expression of the force as a surface integral in Q.) 



A. Conformal coordinate transformations 



The crux of this paper is that the above theory is invariant under conformal coordinate transformations in D- 
dimensional space. These are the angle-preserving transformations r > R, for which the metric transforms as 

9 ' 3 ~* dR* 9km dR> = A (r)5 «' (2 - 13) 

corresponding to local rescaling of distances. The Jacobian determinant of the transformation is thus J = \dR/dr\ = 
X- D (r). 

In a flat (Euclidean) D > 2 dimensional space the group of conformal coordinate transformations comprises the 
rigid transformations (translations, rotations, and reflections) under which the metric does not change; dilatations 
(rescaling) r — > A _1 r, with a constant A, for which g^ = 5ij — * A 2 ^-; and inversions. Under an inversion about a 
sphere of radius a centered at an arbitrary point r a point r is transformed to a point R on the same ray issuing 
from r , with the geometric mean of the distances of r and R from r being a. Explicitly: 
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R = r D + 1 io(r-r ). (2.14) 



The Euclidean metric 5ij then transforms as 



where 



<9i? a 2 

Jy = ~Q^= | r _ ro | 2 (^ - ^iUj). (2.16) 

Here, n is a unit vector in the direction of r — r . The matrix in brackets has eigenvalues 1 (D — 1 degenerate) and 
-1 (non-degenerate). The determinant of Jy, in absolute value, is thus 

T _ a 2D _ \R-r \ 2D 

|r-r |2C - • ( 2 - 17 ) 

All conformal transformations take spheres into spheres (hyperplanes included as spheres of infinite radius). 
Since the scalar potential transforms as ip(r) — > <ys(R) = <p[r(R)], we have V r <p — > Vj?y[r(R)] = (<9r/<9R)V r <,s, from 
which follows that 

(V fl y) 2 = |R " 4 ro|4 (V r y) 2 . (2.18) 

It is customary to use, instead of pure inversions, transformations of the form Pa = IoT(A.)Iq, where T(A) is a 
translation by a vector A, and Iq is the inversion at the origin about a sphere of unit length. These have certain 
advantages: they are connected continuously to unity, and they bring the properties of the conformal group into better 
relief. I prefer, however, to use pure inversions in what follows, as they are easier to handle: they are self inverse, and 
have a simpler transformation Jacobian. 

B. Conformal invariance of the theory 

If r — > R is a conformal coordinate transformation, our action and the field equation are invariant under replacement 
of <p(r) by </j[r(R)], of p(r) by J _1 /j[r(R)], and of the metric <7y(r) by g y [r(R)] (and d D r by d D R in the action). This 
can be checked by direct substitution, but is easier to see as follows: In a general metric space, a theory is conformally 
invariant if its action is invariant under replacement, everywhere in the action, of g^ by C 2 ( r )<?ij (and t hus of g' 



by C~ 2 9 ij j an d of g by C 2D g ), an d of p by C°P [because of the facto r g^ 1 ^ 2 in the definition of p-see eq.(2/7)] It is 
evident from expression ( |2.5D for the action, or from the field equation ^. 8[ ), that ours is indeed a conformal theory by 
this definition. This implies conformal invariance in the above sense, evident by applying first a conformal coordinate 



transformation, under which ip(r) — > y>[r(R)], p(r) — > p[r(R)], and g w (r) — > J 2 ' D g ij [r(R)] [see eq.(2.13)]. The 



action, being a coordinate scalar, is invariant. Now transform the metric back by multiplying it by the conformal 
factor £ 2 = ,J 2 I D , and p by C°- The action remains invariant by virtue of its CI. The net result is that the action is 
invariant under the transformation described at the head of this sub-section. 

It follows from this that if ip(r) solves the field equation for the source p(r) and metric ^(r), then <£(R) = 
y>[r(R)] solves it for the source p(R) = J -1 (r)/j[r(R)], with the same metric g ;j [r(R)]. Clearly, equipotential surfaces 
are transformed into equipotential surfaces. Also, field lines go to field lines, because they are perpendicular to 
equipotential surfaces and angles are preserved in the transformation. Charges are preserved in the transformation; 
i.e., the total charge in a certain volume is the same as the transfo rme d charge in the image of that volume. 

The tracelesness of the stress tensor follows by employing eq.(2.9) with SSf = for 8g tj = e(r)g ijl with e an 
arbitrary, (infinitesimal) function. 

The application of such conformal invariance is standard in the linear, two-dimensional case (in electrostatics, in 
potential-flow problems, etc.). In D dimensions such application have special value because the symmetry is our only 
recourse in solving some of the problems in this strongly non-linear theory, as I do in sections III-V. 

The covariant Laplace (free) action, propto J g 1 / 2 g t3 ip i ip ^ d D r, is not CI in the above sense, but can be made so by 
adding to the above action a term proportional to Rip, with R the scalar curvature, and taking (p to have non-zero 
dimension, so that it transforms as <p — > A^- / 2-1 )^ (see e.g. H). The Euclidean Laplace theory thus becomes CI 
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with <p of dimension D/2 — 1 (as the term with the curvature vanishes), but then the CI of the interaction term 
J g x / 2 ptp is lost. What is special about our theory, and what leads to the applications below, is the fact that it is a 
CI theory in the presence of sources. 

Hereafter I confine myself to the Euclidean case. In curved spaces that are conformally flat, such as maximally 
symmetric spaces, conformal invariance implies the existence of coordinates in which the theory takes the Euclidean 
form, with g i:j replaced by 8ij everywhere. 



C. Asymptotic behavior of the potential 

If the sources p are contained within a finite volume, and sum up to a total charge Q ^ 0, the field becomes radial 
at infinity, and, applying Gauss's theorem to the field equation for a sphere of a large radius, r, we find asymptotically 

V(^«s(QG)|GQ| 1/(D - 1 V- 1 n r . (2.19) 

Here s(x) — sign[x), and n r is an out-pointing, radial unit vector. The potential is then logarithmic for any dimension. 
When Q = 0, the asymptotic behavior of (p is determined by higher multipoles. Typically, a dipole potential dominates 
asymptotically, and has the form ip oc z/r 2 (see below) with z the axis along the dipole. Outside a spherical distribution 
of zero total charge the field vanishes. 



D. Scaling properties 

The field equation enjoys a two-parameter family of scaling invariances: If p(r) solves the equation for a source 
p(r), then, for any two constants a and b, <p(r) = a~ 1 \a/b\ d tp(br) solves it for p(r) = apibv), where d = D/(D — 1). 
When b D = a > 0, so that the total charge remains the same, the scaled potential is <p(r) = ip(br). 

It follows then that the potential, the electric field, the forces, etc. scale simply with charge: If p — > ap, than 
<p — » s(a)|a| 1 ^ D_1 V: an d forces (which scale like gV<p) F — ► |a| d F. These quantities also scale with system size. 



III. EXACT SOLUTIONS FOR THE FIELD 



Onl y fe w charge configurations with exact solutions are known for the general case with an arbitrary form of p(x) 
in cq.( [l.l[ ) 111]]. In particular, there is a closed- form solution for any configuration with one of the D one-dimensional 
symmetries: plane-parallel, cylindrical,..., spherical: By applying Gauss's theorem we get for the present theory 

g (3.!, 

where R is the only coordinate on which ip depends, and s = for the plane-parallel case, s = 1 for the cylindrical 
case, etc.. For a spherical system s = D — 1, and we have 

dp = [Q(r)] 1 /^ 
dr r 

where Q(r) is the accumulated charge at spherical radius r (here and in the rest of the section I use G = 1). In the 
plane-parallel case 

g = [a D S(.)/2] 1 /(^- 1 ), (3.3) 

where S(z) is the total surface density to the left (small-z) of z minus that to its right. 

These solutions, and others, may be used to generate new ones by applying conformal transformations to the 
corresponding charge configuration. Some examples follow. 
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A. Two opposite point charges ±q at ri and V2 



Start with a point charge q > at ri and a spherical shell evenly charged with charge —q, centered at ri, and 
having a very large radius (infinite in the limit). Upon inversion about a sphere of radius a = |ri — r2 1 centered at r 2 
the large spherical shell is transformed into a point charge — q at r 2 , and the charge q stays at rj. The potential for 
the original, spherically symmetric system, is <p(r) = q 1 /^ £)_1 ^n|r — ri inside the spherical shell, and ip = outside. 
It transforms into 

<p(r) = q^/^-^ln^^r, (3.4) 
|r - r 2 | 

(after subtraction of the constant q 1 /( D ~ 1 Hn\v2 — ri|); this applies everywhere. Interestingly, this potential is just the 
sum of the potentials of the two individual charges. This happens to be the case only for two opposite charges. It 
holds neither for two charges that are not opposite, or for more then two charges. 



B. The pure-dipole field 



Asymptotically, at r S> I, where I is the dipole separation, the potential in eq.(3.4) becomes 

^-gVP-^i., (3.5) 

where z is the dipole axis (positive charge to the positive- z side). This is potential for a pure dipole of strength 
qi D ~ x . It describes the field everywhere in the limit £ — > with q^ D ^ constant. (For D > 2, a standard dipole with 
I — > and ql finite does not contribute to the dipole field, due to self-screening effects.) The pure dipole potential is 
a vacuum solution of the field equation that is obtained from another vacuum solution: a constant-gradient field; the 
latter has the potential ip oc z which transforms into z/r 2 . 

The dipole field has a field strength, \Vip\ oc r~ 2 , that depends only on r-not on the angular coordinates. Thi s is a 



well-noticed property of the dipole field in two dimensions. Here it follows directly from the transformation law(2.1S) 
for |V<£>|, and the fact that the dipole field is obtained by inversion from a constant-gradient field. 

For a bounded density distribution of a vanishing total charge the asymptotic behavior of the field is, gcncrically, 
dominated by a dipole field Az/r 2 . I have not been able to express A as a functional of the density distribution. 



C. Point charge in the presence of a grounded sphere 

If, in the above example, we take the charged spherical shell to have a finite radius, we end up with a point charge 
q in the presence of a grounded sphere (as the potential on the original sphere vanishes). By a proper choice of the 
inversion radius and center, the image point charge falls inside, or outside, the image sphere. In the first case the 
potential inside the sphere is that of two opposite charges, and vanishes outside; the tables are turned when the charge 
falls outside the sphere. The charge distribution on the grounded sphere is then easily determined. 



D. Two oppositely charged spheres 

More generally, starting from two oppositely charged (±g) concentric spheres (for which the potential is constant 
in the innermost and in the outermost regions, and is q^'^^Hn r in between) we get the potential field of two 
oppositely charged, equipotential spheres of any sizes, either nested or detached. When the spheres are nested, the 



potentials in the inner, and in the outer part are still constants; in between it is of the form(3.4). When the spheres 



are detached the potential is constant inside the spheres, and is of the form(3.4) outside. 

Starting with two parallel hyperplanes charged with a constant surface density ±S (between which Vy> is constant), 
and inverting about a point halfway between the planes, we obtain two oppositely charged, tangent spheres. The 
potential vanishes inside the spheres; outside we have an exact dipole potential. The charge distribution on each 
sphere-straightforwardly calculated-diverges at the origin and together the charges give a dipole of finite strength. 
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E. Some general comments on potential fields 

Since equipotential surfaces remain so when transformed, and since spheres go to spheres, the equipotential surfaces 
in all the above examples are spheres (all tangent in the case of a point dipole). The field lines are all circles, being 
images of circles or straight lines. For example, for a finite-separation dipole the field lines are all the circles going 
through the two charges. 

There are constraints on the field that can be deduced even when the full field cannot be calculated. As an example 
consider a charge distribution that lies on a circle (with Q = 0). The field Vtp at any point r must be tangent to any 
sphere, of any dimension, containing r and the circle (because the sphere can be transformed into a plane by inversion 
about a point on it). This provides some information on the field of any three-point-charge configuration, or on that 
of a square quadropole. 

Other vacuum solutions of the field equation can be formed by starting from the known, exact solutions of one 
dimensional symmetry (uniformly charged, one-dimensional wire, two-dimensional plane, etc.). As an example take 
a one-dimensional wire in three dimensions with a constant line density a. Working in cylindrical coordinates R, z 
we write the potential as ip = {SuY^R 1 ^ 2 . Inverting about a point off the wire will give the field for certain ring- 
plus-point-charge configurations. Inverting about a point on the wire gives a configuration whose vacuum solution 
is ip cx i? 1 / 2 (i? 2 + z 2 )~ 1//2 . This corresponds to a charge density <r(z) oc z~ 2 (and there appears an infinite opposite 
charge at the origin to compensate the infinite charge of the wire) . 

For a general charge distribution, the field near an arbitrary point, r, away from charges, is conformally related 
to an asymptotic field: Invert about a very small sphere devoid of charges, and centered at r. All the charges are 
transformed into the small sphere, and the r goes to infinity. The asymptotic field of the new configuration is related 
to the field near r in the original configuration. If V<p(r) ^ 0, the image asymptotic field is dominated by a dipole 
term. In an opposite example, look at the field near the mid-point between two equal point charges, where V<p = 0; the 
inverted configuration is a quadropole, with a point charge — 2q flanked by two symmetric charges q; the asymptotic 
field decreases faster than a dipole field. I have not been able to determine this asymptotic behavior. 



IV. MANY POINT CHARGES-GENERAL CONSTRAINTS 



iV-point-charge configurations-comprising bodies whose extent is much smaller than their separations-afford further 
application of the conformal symmetry. Take then a system made of N point charges q\, q N at t\, r N , respectively. 
The information on the dynamics of the system is encapsuled in the energy function E(ri, . . . , r N , q 1 , ...q n ). The energy 
of a charge distribution may be taken as —S. This converges at infinity only when the total charge vanishes. We 
can still use it when Q ^ 0, provided only its changes under change of configuration are needed: only the energies of 
configurations with the same total charge can be compared (see 0). Changes in E may be calculated as changes in 
— S. An infinitesimal change Sp in the charge distribution with no net change in the total charge (J 5p d D r = 0) thus 
produces a change J ip6p d D r in the energy {5p also induces an increment of tp but the field equation implies that 
this does not contribute to the increment of the action). The ./V-point-charge energy, E(ri, r N ), (suppressing the q 
variables) is a special case giving the relative energies of the N point charges in different configurations: we are only 
interested in comparing E values as the charges are moved rigidly to different positions. The energy also diverges 
(logarithmically) near point charges, but these divergences can be subtracted as self energies. I, in fact, treat point 
charges as small but finite bodies, and need not be concerned with such divergences. 

The force on the i'th charge is given by 

F, - (4.1) 

OVi 



I now derive an expression for a virial integral defined in eq.(2.4) that involves only the positions of the charges and 
the net forces on them, but is oblivious to internal forces and structure (the above virial involves integration inside the 
charges). To this end consider the contribution, Vi, of the ith body occupying the small volume v%: Vi = f pr-Vip d D r. 
Write for r within the body, r = r\; + h, where is the center of charge, h is small and will be taken to in the limit. 
We can then write 

Vi = -r i -F i + { ph-Vtp d D h, (4.2) 



where Fj is the net force on the body. The second term in eq.(4.2) does not necessarily vanish in the limit h — > 0, 



because S7<p inside the body diverges in this limit. Let Vy>j(r) be the field produced by the body when it is the only 
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one present. Write V^(r) = V<^j(r) + Vre(r), where k is the increment in the potential due to the presence of the 
other bodies in the system. (Because of the non-linearity k is not just the field produced by all the other charges.) In 
the limit h — > 0, Vifi diverges like | In | 1 (from the scaling properties) but Vk remains finite. Thus, J ph ■ Vk d D h 
vanishes in the limit, and we are left with J ph ■ Wifi d D h. This is just the virial defined for the ith body when it is 



isolated. Thus, from eq.(2.4) 



f ph-Vip d D h^V D {q,) = {dGY^Gq^. (4.3) 
Putting all the above together we finally get an expression for the reduced virial 

V r = - £r« • F« = V D (Q) - V «^) = (dGr\G\ d (\Q\ d 1*1")- (4-4) 

i i i 

Note that the limit of a point charge is gotten from a finite charge distribution p(h) by taking the limit A — > of 
X~ D p(h/ X); this has enabled us to take the limit of J ph ■ Vk d D h to 0. Point-like bodies of finite, higher multipoles 
(dipole, quadropole, etc.) cannot be described in this way (for example, for a point dipole f v , ph dPh remains finite 
in the limit). Our ensuing results are not valid for such bodies. 



A. Scaling behavior of the energy function 



Expression(4.4) implies then that E\ = E(Xr±, Ar n ) satisfies 

f) E 

= Yl x ~ l l Xr < ■ F *( Ar )l = ^IMQ) E M*)l (4-5) 

i i 

Integrating over A between 1 and A we get an important homogeneity property of E(ri, ...,r N ): 

E(Xr u ..., Xr N ) = E(n, r N ) + [V D (Q) - ^ V D (q t )]ln A. (4.6) 

i 

I shall now derive this transformation law of E in a different way, which illuminates better its origin, and which 
will be of further use below. This is based on the invariance of the theory under space dilatations-a fact that also 



underlies the derivation of eq.(4.4) 



Consider first the change in the energy of an arbitrary charge distribution p(r) under a dilatation: p(r) —> p\(v) = 
A~ D p(r/A). In light of the scaling laws described in (IIC), ip(r) — > ip\ (r) = ip(r/X). Integrating SE/S X — j <P\^- d D r 
between 1 and A gives 



E(p x ) - E{p) =lnX p(r)r ■ Vip d D r = V D {Q)ln A, (4.7) 



where Q is the total charge, and I have used eq.(2.4) for the virial integral. This is the result that underlies all our 
findings below. It says that the virial integral is the single system parameter that determines the variations in energy 
under scaling transformations. 

What is the change in E when the charges are moved from positions to positions Ar^? This can be achieved in 
two steps: First apply a space dilatation to the charge distribution. The centers of charge are then moved to the new 
positions. But also, the charges themselves (taken as very small but finite bodies) are dilated by the same factor; this 
is more than we want, as we need to move the charges rigidly to the new positions. After this first step we have from 



eq.(4.7) 



E{Xr u Ar„) = E(r u ...,r N ) + V D {Q)ln A, (4.8) 

where E is the energy of the dilated charges at their new configuration. 

In the second step we dilate each charge separately by the inverse factor to bring the configuration to the desired 
one. The energy change in the second step can be calculated in the limit of very small size for the charges. It is then 
the sum of changes due to separate dilation of the individual charges. This cannot be done when the bodies are not 
much smaller than their separations, and, as before, does not apply if the bodies have finite, higher multipoles. Using 
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eq.(4.7) again for the individual charges yields the a change AE = — In A VoOk). Putting the two together we get 
eq.(4.6). The non-trivial transformation properties of the energy function under — > Ar^ (even when Q = 0) thus 
have to do with the transformation of the self-energies of the point charges. 

Invariance under translations implies that E must be a function of only differences of r,, such that Elji+a) = Efa). 
The derivative of this with respect to a gives 

F = F > = 0. (4.9) 

i 

Similarly, invariance under rotations implies that E depends only on scalars, and that the total moment on the system 
must vanish 

M = J2 r * ® F * ~ F * ® r * = °- ( 4 - 10 ) 



B. Behavior of the energy under inversions 



What does inversion in varian ce tell us about how E changes under inversions, namely, when moving qi rigidly from 
Yi to Ri according to eq.( 2T4 )? We saw above that if the inversion is not to produce a new charge at the center we 
must start with a total charge Q — 0; when Q ^ we may annul it by putting a charge — Q at infinity. Then, from 
conformal invariance, the energy (action) is conserved under a inversion transformation of the charge distribution. 
As in the case of dilatations, such a transformation does not just move the charges to their new positions; it also 
transforms their inner structure; how? When the charges are of very small size their shape change is determined by 
the first derivatives dH/dr, [eq.( 2.16| ) with r and R interchanged]. This describes a reflection about a hyperplane 



perpendicular to n through the body's center, and a dilatation by a factor a 2 /|r< — r n | 2 = R 2 /a 2 , with R, 
As before, bringing the charges back to their original size changes the total energy, which leads to 



iO 



R, 



E(Rx,...,K n ) = E(t 1 ,... ) t ) 



o-£v £ 



(q,)ln 



r?2 

n i0 



(4.11) 



Again, the fact that E is not invariant under inversion of the positions results from the effect on the self-energy of 
the charges. 



The de riva tives of eq.(4.11) with respect to a 2 and to r , at fixed Ri, give sum relations for the forces. The first 
gives eq.(4.4) for the reduced virial, again. The second gives the vector sum relation 



(4.12) 



Equation ( 4.12 ) holds for the forces Fj to which the point charges are subject when at r.;, for all values of a and r . 
Separating the dependence on r and r this equation can be written as 



I - a 2 F + r 2 (l - 2n <g> n) ■ F + 2Vr + 2M-r a = 0, 
where V = ^J r i • F ; - "M<7.)L n = r /|r |, and 

I = rfa - 2 J2( r i ■ F i>i + 2 V Mn- 



(4.13) 
(4.14) 



For eq. fli.15 ) to hold for any r and a we must have separately F = 0, M. = 0, V = 0, and the new sum relation 
1 = 0. The number of such relations totals D + D(D - 1)/2 + 1 + D = (D + l)(D + 2)/2-and tallies with the dimension 
of the conformal group. The first three relations were derived above from transformation properties of E under 
translations, rotations, and dilatations, respectively. Now we see that they all follow solely from the transformation 
properties under inversions. This need not be surprising: locally, all the former transformations can be obtained from 
combinations of inversions. Any set of points in a finite volume can be translated, rotated, or dilated by using a 
succession of inversions alone. 

How do the forces on the charges tran sform under inversion? If Fj and F* are the forces in the old and new 
positions, respectively, then, from eq.(4.11), and choosing the origin at the inversion point (r = 0) 



dE(Hi, ...,Rjv) dr t 



OR 



•F, 



W D (qd% 



(4.15) 
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V. N-POINT-CHARGE CONFIGURATIONS-SOME APPLICATIONS 



Energy functions and forces for certain N point charge configurations can be calculated with the use of the results 
in the previous section. 



A. The two-body system 



The force, F(q 1 ,q 2 ,£), b etw een two point charges q x and q 2 , a distance £ apart [F is positive for attraction) can be 
obtained from expression ( |4.4| ) for the reduced virial: Here Fi = — F2 = F, and the vanishing of the moment implies 
that F points from one charge to the other, so that V r = FI, and we can write 

F = '^[Vo( qi + « a ) ~ V D ( gi ) - V D (q 2 )] = S -^l d -i\G\ d -\\ qi +q 2 \ d - \ Ql \ d - \q 2 \ d ) (5.1) 

[d = D/(D — 1)]. Th is result was derived in ^ for the three-dimensional case in a roundabout manner. (For 
D = 2 expression ( |5.l| ) reduces to the standard two-dimensional linear- medium result F = Gq 1 q 2 /£.) For example, 
for two equal charges (fr = q 2 = q, F = 2s(G)£~ 1 d~ 1 \G\ d ~ 1 \q\ d (2 d ~ 1 — 1). For oppos ite charges: q t = —q 2 = q, 
F = — 2s(G)£~ 1 d~ 1 1 C I rf 1 \q\ d , also to be gotten from the force-transformation law (4.15) starting with one charge at 
infinity, and hence F* = 0. 

The energy function in the two-body case is 

E(r 1 ,r 2 ,qi,q2) = /3i2H r i - r 2 |, (5.2) 

where 

fly ee V D ( qi + qj ) - V fl («i) - V D ( qj ). (5.3) 



B. The three-body system with vanishing total charge 



Consider three charges q { at r i} with qi + q 2 + q$ = 0. Perform an inversion with one of the r^s as center, say 
r3. Then, 53 is transformed to infinity, and q\ t 2 are transfo rmed to Ri^- The force on qi, say, can be calculated in 
the new configuration from the two-body force formula eq.(5.1). From t his t he force Fi in the original three-body 
configuration is calculated by employing the force-transformation law eq.( 4.15| ) to obtain 



Pi = A 



i"2 - ri 



■0i 



r.3 - ri 
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(5.4) 



(This could also be derived from the above constraints on the forces F = 0, M — 0, V = 0, 1 = 0.) Interestingly, 
the force is the sum of the two forc es that would have been exerted by qi and (73 separately, the non-linearity 
notwithstanding. Integrating eq.( |5.4|) over ri we get the explicit form of the three-body energy function (for the 
zero-total-charge case): 



E(r 1 ,r 2) r 3 ) = ln[r^ 



Pl2 r 013 P23} 

12 '13 ' 23 J 



(5.5) 



C. The virial theorem 



I now derive the analogue of the standard virial theorem, (Ek) = —{E p }/2, relating the mean kinetic energy and 
mean potential energy of an iV-body system held together by Newtonian gravity in three dimensions. Consider a 
bound system made of any number of point charges q t , of masses m il moving under the sole influence of the ip field 
they produce (other forces may act inside each body to hold it together). The center-of-mass acceleration of each 
body t\ =F i /m i , with the Fj satisfying eq.fl4.4p. Now 
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I i i i 

The first term vanishes in the stationary case (or its long-time average vanishes for a general bound system); the 
second term equals twice the kinetic energy Ek — M(V 2 )/2, where M is the total mass of the system, and (V 2 ) is 
the mean square velocity. Together with eq.( [4.4[ ) this finally gives the desired virial theorem 

2E k = M{V 2 } = V D (Q) - - (dGr\G\ d (\Q\ d - £ |g,| d ), (5.7) 

i i 

by which the mean-square velocity depends only on the charges. This is exact for a stationary system; for a general 
bound system the long-time average of the left-hand side has to be taken. 

This relation has been used to estimate the mass of large-scale, cosmological, galaxy filaments, which are approxi- 
mately two-dimensional Newtonian systems |(|. It has also been used to determine the masses of roundish galaxies, 
in the modified dynamics, from their observed velocity dispersions, when typical accelerations in them are very law 
(see Q and references therein). In these cases, where g 4 are the constituent masses, and where N ^> 1, ^ | | rf 
can be neglected as it is smaller than \Q\ d — M d by a factor ~ jy- 1 / 1 , - 1 ) . In the limit N — > oo 

(V 2 ) = ^l(GM) 1 /^-!). (5.8) 



D. Symmetric configurations 



The force on bodies in some symmetric configurations can be calculated from expression (4.4): Consider a config- 
uration comprising a charge go at the center, and n equal charges, g, at positions r.; that are equivalent with respect 
to the center-equivalent in the sense that each of the points, r^, can be interchanged with any other by an element of 
the symmetry group of the system, H: a rotation, a reflection about the center, a reflection about some hyperplane 
through the center, or combinations thereof, that is also a symmetry of the system. Examples are the corners of 
a rectangular hyper-box, the vertices of any perfect solid of dimension D or less (such a s a perfect polygon, any 
hyper-cube, etc.), the vertices of polygonal prisms of different types, etc.. Clearly, the g-charges are then all at the 
same distance from the center, call it r. Also, they are subject to equivalent forces F^; to wit, each of the F^s can be 
transformed to any other by an element of the symmetry group. In particular, the radial components of these forces 
(F • r)i are all equal, because scalars are invariant under the point group. (The force on the charge at the center 



vanishes.) This common value can be deduced from eq. ( 4.4 ) (the symmetry automatically insures that F = 0, M = 0, 
and I = 0): 

F • r = -~[V D (q + nq) - V D (q )] + V D {q). (5.9) 
n 

When the forces are radial-as when the point are the vertices of a perfect solid-the full force is obtained since in this 
case F ■ r = —Fr (F is positive when F acts towards the center). 



In the limit n — > oo, nq — > Q, eq.(5.9) gives the force on the elements of a spherical shell of any dimension smaller 
than D (e.g. a ring) having radius r, total charge Q, and a charge go a t its center. The force-per-unit-charge on the 
shell is 

F=-^{V D (q Q + Q)-V D (q )}. (5.10) 
Qr 

Additional results are described in appendix A. 



E. Point charges in the presence of conducting boundaries 



When (equipotential) bodies of infinite conductivity are present, the full conformal symmetry enjoyed by the N- 
charge problem is destroyed (the boundaries remain in place when applying the transformation to the charges). So, 
for example, homogeneity is always lost, and the total force now does not vanish in general; rotational symmetry 
sbout an arbitrary center is lost, an so the total moment does not vanish; etc.. Some symmetry may however be 
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left, in which case the corresponding identities are still valid. If the arrangement of conductors is invariant under 
translations in a certain direction, the component of the total force on the charges in this direction vanishes. If the 
conductors are spherically symmetric about some center, the total moment on the charges with respect to this center 
vanishes, and so on. An interesting example involves boundaries that are invariant to rescaling about a certain point 
(taken at the origin). This happens when for every point, r, on the boundary, Ar is also on the boundary, for every 
A > 0-the boundary is an arbitrary-cross-section cone (including a hyperplane, a corner, etc.). In this case, expression 



(4.4) for the reduced virial holds with respect to the origin (but 1 = does not); Q now is the total charge including 
that on the conductors. When the conductors are grounded they automatically take up a charge that makes Q = 0. 
Thus, for example, the force on a single charge q, in the presence of an arbitrary, grounded, conic conductor is always 
subject to a force F satisfying 

F-r = V (q). (5.11) 

This can be extended to conductors in the shapes of discs or spherical caps, as they can be transformed into half-planes 
by inversions. 



VI. OTHER CONFORMAL ACTIONS? 

Rotation invariance dictates that an action containing only first derivatives of ip is a function of (Vp) 2 , and scale 
invariance further requires that it be of the form S p oc /[(Vt/j) 2 ] p d D r. By themselves, all such actions are invariant 
under the scale transformation p(r) —> cp(r) — X a ip(Xr), with a — (D/2p) — 1, namely, with ip having conformal 
dimension a. If it is to be fully CI, S p must also be invariant under inversion at the origin: ip(r) —> <p(r) = 
(a i 'r) 2a p{a 2 r / V 2 ) . This can be shown not to be the case unless either p = 1 (the linear case), or p — D/2, which is 
our action Sf. To see this note that starting with p oc j- - ^ -2 *')' ( 2p_1 ), which is the only spherical vacuum solution 
of the theory (beside p = const), the corresponding sp is not a solution, unless p = 1, D/2. Similarly, if we start with 
ip oc z, which is a vacuum solution for all S p , the corresponding p is not, unless p = 1, D/2. 

Field theory lore has it that scale-invariant theories tend to be CI, but this is not a theorem (see e.g. fH| , and 
references therein). The above nonlinear theories constitute counter-examples. 

For the interaction action to be invariant we need p to transform according to p(r) — > yj[R(r)], since p d D r is 
invariant; i.e., p is then of conformal dimension zero. As explained in section IIB, the Poisson action in D > 2 is not 
CI because for the free action to be invariant <p has to have dimension D/2 — 1. We are thus left with S D as the only 
CI action containing only first derivatives in the field part. 

Consider now actions with a field part containing higher derivatives written for curved space 

S=-j g l/2 L{g^,Vp,...,V k p)d D r- J g x l 2 ppd D r. (6.1) 

The field Lagrangian L depends on covariant derivatives of tp: p ;il .... im , 1 < m < k (collectively designated V m p), 
and is a coordinate scalar. Scale invariance alone is tantamount to the homogeneity requirement 

L(gV,XDp, X k V k p) = X D L{g^,Vp, ...,V k p). (6.2) 

This is because coordinate invariance dictates that every derivative in L is contracted with another by one g ZJ . As ip 
has zero dimension, under dilatations, r — > A r, an mth covariant derivative is multiplied by X m , which is the same 
as multiplying g ij by A 2 , and conformal invariance tells us that L(X 2 g ij , ...) = X D L(g ij , ...). 

Consider, hereafter, the Euclidean version of the action (where <? y is put to and T> to d). As stated above, scale 
invariance does not insure full CI. The quadratic Lagrangians-leading to linear field equations-with p of zero scaling 
dimension: L — p>A D ^ 2 p (for even D) are CI. The non-linear actions with p of zero scaling dimensions are probably 
not. In fact, I have not been able to find any that is ( without having a general proof that none is). For example, I 
have shown that all the Lagrangians, in even D dimensions, of the form 

L = [(vv) 2 r(A^) fe , (6.3) 

with m = D/2 — k, and with k = 1 for D > 4; k = 2 for D > 4, or k > 3 for D > 2k, which are scale invariant, are 
not CI. 

The homogeneity condition( |6.2j ) implies that the field equation necessarily has a vacuum, spherically symmetric 
solution of the form p = Aln r. The Euler-Lagrange equation can be written in the form diJt — p, where Ji is a ve ctor 
that is a function of the first 2k — 1 derivatives of p, with the homogeneity property [easily derived from eq.(|6.2D] 
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JiiXdip, X d ip) = X Ji(dip, (6.4) 



When i/3 depends only on the radial coordinate r, J has only an r component, and, from eq.(6.4), for ip = Aln r, 
J r = r~( D_1 ) j(A). Since in the spherical case diJi = 7 ,_ (- D_1 '9 r [r- D ~ 1 J r ], clearly ip — Aln r is a solution. The 
coefficient A is determined via the Gauss theorem: j(A) oc Q, where Q is the total charge at the center. If L is 
homogeneous in ip of degree /?, then J is homogeneous of degree f3 — 1 in and then the coefficient A is given by 
Aocsi 5 n(Q)|Q| 1 /(/3-i). 

The purely logarithmic potential outside a single spherical body is, however, valid for only very specific density 
runs. When the Lagrangian depends on derivatives up to the fcth, the generic, spherically symmetric, vacuum solution 
is characterized by 2k constants, which are determined by the exact density run in the central body. In general, the 
potential diverges at large r as a power in r. The notion of a point charge is thus not useful as the external solution 
does not depend only on the total charge. We cannot even speak of "the field of a point charge" as this is not well 
defined. For example, the quadratic theory with if of zero dimensions, with L oc tpA D / 2 ip + Apip, with the field 
equation A D ' 2 ip oc p, has D independent spherically symmetric vacuum solutions of the form <p = const., In r, r ±Q , 
with a = 2,4, , ...,D - 2. 

VII. MULTI-POTENTIAL THEORIES 

The above theory is straightforwardly extended to describe K (coupled) scalar potentials, ip a , which couple to K 
types of charges with densities p a (a = 1, K). The action is 

with a a = (Vipa) 2 , and G = 1. Conformal invariance is now equivalent to homogeneity of C: 

C(Xa a ) = X D/2 C(a a ). (7.2) 

(One could actually generalize further by taking a a b — Vip a ■ V<Pb as variables, but I keep to the simpler form.) The 
K field equations are 

V- \p a (ai,... )V(p a ] =a D p a , (7.3) 

with ji a = dC/da a , I < a < K . The spherical vacuum solution for a system with total charges q a is 

fa = siqJQW-Vln r, (7.4) 

with Q a > 0, giving a a = (Q a ) 2 ^ D_1 V -2 ; s(q a )Q a may be viewed as the asymptotically observed charges. They 
are determined from the actual charges, q a — J p a , as follows: Inserting the above form of ip a in the Gauss theorem 
p, a (cti, ...)d<p a /dr = (7 Q r _ ^ D_1 \ and making use of the homogeneity property of the [i a (derived from that of £): 
p(Xxi, Xx2, ■■■) = X^ ^ 2 ^ 1 ^ /i a (xi, X2, ■■■), we get the K equations: 

^[Ql /{D - 1 \Ql /{D - 1 \..]Ql /{D - 1) = ka\. (7.5) 

All our results in the previous sections can be carried through mutatis mutandis to this, more general, case. The 
virial integral, which controls the change in system energy under dilatations, can now be shown to be given by 

V. / $>r. V^ a ^=^£[Q 1 ^- 1 \Q^,...,Q^~ 1 ']. (7.6) 



Using the homogeneity of £, eq.(7.2), by which C{x\, xk) = (2/-D) J2a x a^a), and with eq.(7.5), we can also write 



V=^£k^ /(D ~ 1) - (7.7) 

a 

This expression for V is to replace V D (q) in all our results (remember that the g a s and Q a s are charges of different 
types of t he same body). For example, the powers (3ij appearing in the two- and three-body energy functions defined 



in eq.(5.3) are now given by 
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fki = —y-{£[Q 2 i /(D1) > ■■•] - C[{Q\f'^\ ...] - Cmifl^, ...}}, (7.8) 

where Q l a a = 1,...,K are the "asymptot ic" charges for the point body i, and Q a are those for the two bodies i,j 
taken together; i.e., as calculated from eq.( |7.5| ) with the charges q a = q l a + q^. 

As a special case, assume that the theory is invariant under rotations in the internal space of tps, i.e., under 
fa —> fia = Oabfb where O is an orthogonal matrix (the charges are then rotated by the same matrix, leaving 
^2aPaVa invariant). This means that C must be a function of X) a a ai ana - ^ ne required homogeneity then dictates 
that 




(the constant in front is chosen to match the single-field case) . 

Equation ((7I]) can now be solved to give the asymptotic charges as 

Qa = \q a \ D - 1 q- {D - 2 \ (7.10) 

where q = i^2, a q 2 ) 12 is the root- mean-square over all the charge types of a body. The virial now takes the single-field 
form, only with the system's single charge replaced by its root-mean-square over all the charge types: V = V D {q). 



VIII. VECTOR AND HIGHER-FORM THEORIES 



Maxwell's electromagnetism is governed by the action 

S = ~J 9 1/2 F^F^ d D r + J g 1/2 J f "A f , d D r. (8.1) 

It describes the electromagnetic field F^ v = A v ^ — A^ M , derived from the vector potential A^, in the presence of 
conserved currents J^. The theory is gauge invariant, and is conformally invariant in four dimensions only (i^" = 
F a pg a,f 'g f)l ' , and J M = <7 _1 ' 2 i /i , where the vector density contains only matter degrees of freedom, but not the 
metric). In a vein similar to that in our treatment of the scalar case we can construct non-linear, CI generalizations 
for D > 4: Take the field action to be 

S f = - J g- 1/2 C(F, lv ,g, lv ) d D r, (8.2) 

where £ is homogeneous of degree D/2 in F: L{j]F ilVl g llv ) — i] D/2 C(F^,g^) (and, of course, is a coordinate scalar). 
It is CI because multiplying g^ v by A(r) in C is tantamount to multiplying F^ v by A _1 (r) (every two lower-case indices 
in F^ must be contracted using one g^)- A factor A~ D/2 is then pulled out of C to cancel the factor from g 1 ^ 2 . In 
D > 4-dimensions we may take 

C oc {F^F^) D '\ (8.3) 

which gives a CI, vector theory when coupled to the currents as above, but there are others. For example, in eight 
dimensions 

C cx F^F^F^F^ (8.4) 

is such a theory. 

More generally, we have linear, gauge-invariant, CI theories in even D dimensions involving an antisymmetric- 
tensor gauge potential of rank n = D/2 — 1 (an n-form potential), 4 aii .., ian The field tensor -ff ai ,...,a„ +1 is the 
totally antisymmetrized derivative of A, and the current density is also an n-form, in analogy with the Maxwellian 
case. The field Lagrangian is C oc -ff Ql ,...,Q„ +1 -ff Ql ' " ,Q " +1 , and the interaction Lagrangian is A ait .„ iCtn J a i. I n 
dimensions higher than 2(n + 1) we get a CI theory by taking C{H) that is homogeneous of degree D/(n + 1) in H, 
e.g. C oc (i?Q 1 ,...,Q, l+1 -ff Ql ' - ' a " +1 ) £ ' /2( " +1) - There is no known linear, CI generalization in D > 2(n + l)(see e.g. g). 

Specialize now to flat spaces, I find that, as in the scalar case, the effect of rescaling on the energy of some current 
distribution is given by 
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E x - E = V In A, 



(8.5) 



where E\ is the energy of the rescaled current distribution. Again, V can be written as a surface integral. For example, 
for the vector-potential case, where the rescaled current density is A"'"^ 1 ' J M (r/A), we have 

/BC 
da„ ( r »£-2rPA a . P7 — ). (8.6) 

This can be shown to vanish for configurations with spatially bounded currents (bounded in space and time in higher- 
D Minkowski spaces). We do not have the appropriate analogues of point charges in the scalar case, with finite V, 
for which to calculate and manipulate TV-point energies. The situation is akin to having, in the scalar case, a system 
of point bodies of null charge but a finite higher multipole. As explained in section IV, our treatment of point-charge 
systems does not carry to such systems. 

We can still employ the CI to create new solutions of the non-linear theory from other known solutions. For example, 
the vector potential = (l/2)B afl r a , with B all constant and antisymmetric, gives a constant field = B^, and 
is thus a vacuum solution of all the above vector theories in D dimensions. Inversion at the origin, under which 
Afj,(r) — > r~ 2 A fJi (r/r 2 ) (A^ is of dimension one in these theories) gives A^ = (l/2)B afi r a r , which must them also 
be a solution. 



IX. DISCUSSION 



It was pointed out to me by David Kutasov (private communication) that our results for the classical field theory 
evoke, and might well be related to, results known to hold in conformal field theory. Perhaps there is such a conformal 
QFT whose classical limit is our theory. Comparing with the definition of a conformal QFT as given e.g. in 0, to 
make such a connection we will have to identify the so-termed "quasi-primary" fields of the QFT with e lqip ^ r \ where 
ip(r) is now a quantum field. The TV-point-charge energy function is the classical limit of the correlation function of 
these quasi-primary fields 

E(n, iv ) = ln[{e iqi ^ ■ ... ■ e ^f{r N )^ (g 



and the constraints I found such as (4.4)( 4~T^ ) are the Ward identities for the correlator. What I then proved amounts 



to showing that the (anomalous) dimension of the operator e lqv ^ is V D (q). 

In the vector- and higher- form-potential case we cannot identify such an in finite set of "quasi-primary" fields. 
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A. SOME MORE EXAMPLES OF FORCE CALCULATIONS 



The forces on the charges can also be calculated for symmetric configurations of the following type: Start with the 
symmetric configuration described in section VD, with an even number of sites, n. Put the charge at the center to 
0. Place an equal number of positive and negative charges ±g at the sites in such a way that all charges of the same 
sign are equivalent: Every two charges of the same sign can be interchanged by a point symmetry that does not mix 
charges of a different sign. The corners of a 3 — D rectangular box, for example, can be decorated in three inequivalent 
ways that satisfy the above. The vertices of a perfect 4m-polygon can be decorated in two ways: alternating charges, 
and in a two-pluses-two-minuses pattern. Here again the forces on all charges, positive as negative, are equivalent 
(can be transformed to each othe r by elements of the full, site point group H). Again, F • r takes the same value for 



all the charges. Using expression (4.4) we get 

F-r = V D {q). (A.l) 

A charge configuration of the above description results when we apply the method of images to the problem of a 
charge q in a region bounded by two intersecting, grounded, D — 1-dimensional hyperplanes. When the angle between 
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the hyperplanes is a — ir/m, the field in the region bounded by the planes is the same as that of a system of images. 
This has m pairs of charges iq arranged alternately on the vertices of a polygon in a configuration as above (the 
polygon is not perfect bu t has edges of alternating lengths). Equation (A.l) thus gives the radial force on the charge; 
it is a special case of eg. ( 5. 11 ). 

Now an example of the use of the constraint 1 = 0, where I is defined by eq.(4.14). Consider a planar, perfect 
polygon of n equal charges q (a uniform ring in the limit n — > oo), and a charge —nq on the symmetry axis of the 
polygon, a distance £ from the origin at the polygon's center. We want the force F on the large charge (which acts 
along the axis), and the force f on the small charges. From I = we have 



F 



2V D (ng)£ 
r 2 +£ 2 ' 



where r is the radius of the polygon. The axial component of the force on the small charges is 
component, f r is gotten from V = 0: 



(A.2) 

-F/n, and the radial 



fr = -r~ 



V D {q) + -V D (nq)- 2 



f 2 



( 2 



(A.3) 



(F and f r are positive when towards the center). Since this configuration is obtained by inversion from one in which 
the large ch arge is at the center of the polygon, the above results also follow by applying the force-transformation 
formula eq. ( 4.15| ) to the results of the symmetric case. Inversion about a point placed on a uniformly charged ring 
(limit of a polygon) transforms it into a line with charge density p{x) — po[l + (^A^) 2 ] -1 ; an d the point charge at the 
center can be moved to a point at an arbitrary distance in the symmetry plane of the wire. 
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